Mass Shift and Width Broadening of J/ ip in hot gluonic plasma from QCD Sum Rules 
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We investigate possible mass shift and width broadening of J/tp in hot gluonic matter using QCD sum rule. 
Input values of gluon condensates at finite temperature are extracted from lattice QCD data for the energy density 
and pressure. Although stability of the moment ratio is achieved only up to T/T c ~ 1.05, the gluon condensates 
cause a decrease of the moment ratio, which results in change of spectral properties. Using the Breit-Wigner 
form for the phenomenological side, we find that mass shift of J/ip just above T c can reach maximally 200 
MeV and width can broaden to dozens of MeV. 
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Heavy quarkonia have been regarded as a very useful probe 
of quark-gluon plasma (QGP) which can be created in rela- 
tivistic heavy ion collisions. Since Matsui and Satz argued 
HI] that the J 'ftp suppression could be a signature of QGP 
formation in heavy ion collisions, extensive works have been 
performed on the subject both experimentally ||2|] and theo- 
retically J3|. However, in contrast to early expectations, re- 
cent lattice calculations suggest J/ ip can survive up to at least 
T ~ 1.5T C JH,§,|^]. Hence, change of spectral properties, 
which cannot be seen within current resolution of lattice cal- 
culations, may exist in QGP and can reflect the properties of 
the strongly coupled QGP, at temperature not so higher than 
T c . One should note, however, that in such a temperature re- 
gion, the system is highly non-perturbative one |7|], and any 
analyses should consistently treat the non-perturbative aspects 
of QCD. In this respect, QCD sum rule is a suitable approach 
that can be used for analyses of hadron properties at these 
temperatures. 

In the present paper, we investigate the behavior of J/ ip in 
gluonic plasma slightly above T c using QCD sum rules |8]. 
The QCD sum rule has been extensively used for studying 
in-medium properties of both light and heavy hadrons as a 
reliable and well-establish method. For heavy quark systems, 
it is more reliable because all relevant condensates are known 
and their temperature dependence are easily extracted from 
lattice QCD. Here we follow the method used in analyses of 
J/ip in vacuum [9] and in nuclear matter Hoi [Till and study 
J I tp at finite temperature within the quenched approximation. 

We start with the time ordered current-current correlation 
function 

(q^-q 2 g^)IL(q)=i [ d 4 x e^ x (T[f{x)f (0)]) T (D 



where we take j p — cj^c for J/ip. In this work, we set both 
the medium and cc at rest so that q = and tl(q 2 ) becomes 
the longitudinal part of the polarization tensor. For a large 
Q 2 = —q 2 = —oj 2 ^S> 0, the correlation function can be ex- 
pressed through OPE as H(q 2 ) — C n (O n )T with C n and 
O n being the perturbative Wilson coefficients and operators 
of mass dimension n, respectively. In heavy quark systems 
such as J /ip, the expansion can be entirely expressed by glu- 
onic operators J§,@, T^, HI- Here we assume that the relevant 



information is contained in the local operators. Note that, in 
the deconfined phase, this may have to be remedied since non- 
local contributions coming from the non-vanishing Polyakov 
loop might become important [12]. The expectation value is 
taken at finite T > T c and the temperature dependence can be 
imposed only to the expectation value of the gluonic operator 



in the case of m c > T and T 2 < Q 2 tlj. In previous QCD 



sum rule works, Hashimoto et al. calculated the Wilson coeffi- 
cient of the scalar gluon operator at finite temperature I14I1 . In 
the works of Furnstahl et al. 11511 . scattering contribution was 
calculated. However, none of these works included the con- 
tribution from the gluon operator with spin nor the changes 
of the gluon operators systematically extracted from the lat- 
tice calculations. In contrast, we have included all the lowest 
non-vanishing operators, have extracted the full temperature 
dependence of the gluon operators from the quenched lattice 
data, and have modeled the phenomenological side consistent 
with the quenched assumption. Hence, this work marks the 
first systematic application of QCD sum rules slightly above 
%. 

Following Ref. [9], let us consider the rt-th moment of the 
correlation function (Q~|i, 
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From the OPE side, up to dimension four, this moment can be 
expressed as 

M n {Q 2 ) = + a„(£H + b n {£)<k> + c„(£)^c] (3) 

where £ = Q 2 /Am 2 is a dimensionless scale factor, A n , cz n , 
b n and c n are the Wilson coefficients corresponding to bare 
loop diagrams, perturbative radiative correction, scalar gluon 
condensate and twist-2 gluon operator, respectively fioh . Note 
that there appears the additional twist-2 contribution when we 
consider the medium expectation value. The Wilson coeffi- 
cients are listed in Ref. J^l The explicit forms of cp\, 
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to G and cp c = 
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and C are vb - 9(4ra i )4 uo a..u Vc - 3{4m 2yi 
Go = {^G a^G 1 /^) T , and G2 is the twist-2 condensate con- 



tribution defined by (^G% p G v ap ) T = 



3<r ) G 2 

with being the 4-velocity of the medium. These conden- 
sates can be determined from lattice QCD data as follows. 
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The scalar condensate is related to the energy-momentum 
tensor through the trace anomaly. If we take 1-loop expres- 
sion for the beta function of pure SU(3) theory, we get Go = 
G^ ac -^-(e-3p) where G£ ac = (0.35GeV) 4 is the value of the 
gluon condensate in vacuum IU6I1 and the second term comes 
from the trace anomaly lfl7ll . e and p are the energy den- 
sity and pressure, respectively. On the other hand, the twist- 
2 part can be simply related to the energy-momentum tensor 
of the pure gauge theory as T afJ = — G" A G^ A {a ^ (3). 
Hence, recalling that T a/3 = (e + p)u a u^ — pg a/3 , we obtain 
(e + p). The thermodynamic quantities and the 



G 2 = 



effective coupling constant ajT) are taken from quenched 
lattice QCD calculation 111 8l Il9ll . Accounting for the am- 
biguities of a s (T) in the non-perturbative regime, we adopt 
two of results in Ref. lfl9ll for the temperature dependent cou- 
pling constant; one is determined from the short-distant force 
and the other is from the screening part of the large distant 
part. The former does not depend on temperature at very 
short distance and takes its maximum values at some distance 
^screen which decreases with increasing temperature. We use 
the value at this distance and denote it as a qq {T) following 
Ref. 01911 . To obtain the temperature dependence, we fit the 
lattice data point (Fig. 6 top in Ref. IU9IP by Bezier interpo- 
lation, which results in a q „ (T r ) = 0.626. The latter, which 
we denote a(T) as in Ref. 11911 . has a similar value, but error- 
bars are still too large especially near T c . We use the two-loop 
expression of the running coupling constant with a set of pa- 
rameters given in Ref. I19L This gives a(T c ) — 0.47. The 
extracted gluon condensates Go and G2 above but near T c are 
shown in Fig. Q] One thing to note is that, Go decreases to 
less than half of its vacuum value but remains positive near 
T c 112011 . We can see that the tensor condensates have non- 
negligible values near T c . 

The n-th moment in Eq. <(2j can also be expressed as 
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where ph(s) — — tanh 1 



^y) Imn(s) is the phenomenologi 

only 
Gil 



cal spectral function which in general includes not only the 
pole term but also the continuum and scattering part I15ll2lll . 



The scattering term, which also appears in the OPE side and 
contributes with a delta function at zero frequency, could be 
important in the presence of thermal fermion. However, since 
we are considering the gluonic medium and have extracted the 
condensates from the pure gauge theory, we can consistently 
assume that there is no (anti-)quarks which can scatter with 
the current. 

We can put tanh (^y) =1 due to much larger pole mass 
and continuum threshold than temperature considered here. 
Then, hadronic properties such as mass and width are related 
to the OPE side [Eq. ©J by putting a phenomenological func- 
tional form in ImII(s) of the above equation. Here we employ 
a simple Breit-Wigner form 
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FIG. 1: Gluon condensates as a function of T/T c . 



to take finite width into account. Since we are interested in 
the lowest lying resonance of the vector channel, we should 
choose an appropriate order n so that the moment [Eq. 
contains information only on the pole term of the spectral 
function. Following Refs. |@, [loll , we take the ratio of 
the moment r n — M n -\/M n and choose moderately large n 
such that the contribution from the excited states and contin- 
uum can be neglected. Therefore, this ratio should approach 
a constant value at sufficiently large n. However, when n is 
large, contribution from higher dimension operators becomes 
important. At the n value where r n is minimum, pole dom- 
inance and truncation of the OPE are valid and the ratio is 
close to the real asymptotic value, as have been extensively 
investigated in the vacuum sum rule for J/ip |9J]. In this work 
we only consider temperature range in which the same crite- 
rion can still be applied and take the minimum value for r„ 
to be its asymptotic value. Hence, in the practical calculation 
below, we firstly evaluate the appropriate n for various tem- 
peratures by calculating r n |ope- Then, we look for pairs of m 
and r which satisfy the sum rule relation, r„|opE = ?"n|phen.- 
We employ the vegas monte-carlo integration 12211 to treat a 
very sharp peak in the dispersion integral [Eq. (0|] of the phe- 
nomenological side. The relative error in the numerical inte- 
gration is found to be on the order of 10~ 6 for m — 3 GeV 
and r = 1 MeV. This numerical accuracy becomes better as 
r increases, as naively expected. The normalization scale £ 
is chosen as £ = 1. We checked that our result does not 
strongly depend on the choice of £ by varying £ from to 3 
lE3tl . Other parameters of the theory are taken from Ref. Iioll . 
a s (8m 2 ) = 0.21 and m c = 1.24 GeV. We did not do fine 
tuning of these parameters to adjust the vacuum mass of J/ ip 
since our interest is in the change of mass and width induced 
by the hot medium. Figure [2] displays the ratio of the mo- 
ment calculated from Eq. ([3]) using a(T) for G2. We can 
see that the ratio have a stable point from vacuum to 1.05T C 
but the stability is no longer achieved beyond 1.06T C . Also 
as seen from the figure, the stable point shifts to larger n as 
temperature increases. However, Eq. OJ becomes worse as 
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FIG. 2: Ratio of the moment from the OPE for various T/T c . 
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FIG. 3: The moment ratio of the phenomenological side as a function 
of F for various masses. Left figure shows n = 9 case and right one 
shows n — 14 case. 



n increases, because the Wilson coefficients increases with n 
0. This can be improved by increasing £, but the stability 
holds only up to 1.06T C even for £ = 3.0. If we use a qq (T) 
for G2, the stability becomes worse due to its larger value at 
this temperature region. In the £ = 1 case, there is no sta- 
ble point for T = 1.05T C with a qq {T). This lack of stability 
does not necessarily mean dissociation of J/ip but shows a 
breakdown of our approximation. The reasons for the break- 
down are twofold. One is the lack of convergence of the OPE 
in Eq. (01. This can be improved by including higher dimen- 
sional operators [24]. The other is physical one. Since the 
non-perturbative part largely decreases above T c , perturbative 
contribution will become more important. Hence, in order to 
study higher temperature region, we will need to improve the 
phenomenological side to be more consistent with the OPE 
side, which can be accomplished by a temperature dependent 
continuum contribution. Then it will lead to n-independent 
results for physical parameters until the J/tp really dissolves. 

Before going to results of m and T, it is useful to see a 
feature of the phenomenological side. We depict the moment 
ratio of the phenomenological side based on the dispersion in- 
tegral (|4]i in Fig. [3] We choose two n values which correspond 
to the stable points at vacuum and T = 1.05T C , respectively. 
We see that the moment depends on T very weakly. For the 
larger n, the dependence becomes slightly stronger. On the 
other hand, the moment of the OPE side (Fig.O shows about 
2 GeV 2 decrease from vacuum to 1.05T C . Hence, the width 
must become very large to achieve 2 GeV 2 reduction of the 
moment ratio if change of the mass is small. The process to 
determine the mass and the width is nothing but evaluating the 
intersection between the stable points in Fig. |2] and the curves 
in Fig. [3] Because of the monotonic behavior of the r n as a 
function of m and T, mass satisfying the sum rule takes its 
minimum value in the r — * limit, i.e., magnitude of the 
mass shift becomes the largest if width stays constant. In this 
limit, mass is given by a simple relation m 2 = r n |oPE — 4ra 2 £ 
J9I]. We cannot determine both mass and width only with the 
sum rule because it provides only one equation with respect to 



two unknown quantities, m and T. The situation is similar to 
light vector meson l25ll . However, we can extract a relation be- 
tween the mass shift and width by fixing the mass firstly and 
then solving r„|opE = fnlphen. f° r T because the monotonic 
behavior of the r n guarantees the unique solution. 
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FIG. 4: Relation between mass shift and width at finite temperature. 
Left figure shows the case of a qq (T) and right one shows the case of 
a(T). 

In Fig. 2J we display the result of relation between mass 
shift 5m — m me dium — m vacuum an d width which satisfy the 
sum rule. The result gives a clear, almost linear relation be- 
tween r and 5m. The difference in a s (T) appears as 10-50 
MeV difference of the mass shift at T = and about 10 MeV 
difference of the width at 5m = 0. 

Finally we plot the two extreme cases, mj^ for T — > and 
r for 5m — > 0, as a function of temperature in Fig. [5] In both 
cases, the change is almost linear with temperature. We stress 
that, however, these results are extreme cases. Our results 
show there must be notable change of mass or width, or both 
of them. Hence, once either mass or width is estimated by 
other methods, one can obtain the other through the relation 
given in Fig. @] For example, a pQCD calculation can give 
thermal width 0260. This does not show large (~ 100 MeV) 
thermal width near T c . Then the large mass shift is expected. 
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FIG. 5: Temperature dependence of the mj« (left vertical axis) in 
r — » limit (max. mass shift) and Tj/^ (right vertical axis) in 
5m — > limit (no mass shift). Mass and width are indicated by open 
symbols and closed ones, respectively. 



Such mass shift is in fact expected, as the sudden reduction 
of the asymptotic value of the potential just above T c seen in 
a lattice QCD 127I1 will inevitably lead to lowering of bound 
state energy on that potential 1 28 , 2^] . Although the accuracy 
of the lattice MEM method is not enough for clear compar- 
ison, the spectral function calculated from a potential model 
motivated by a full-lattice QCD shows a large shift of the peak 
Q. A recent full lattice QCD also shows the shift of the J/ip 
peakjl. 

In heavy ion experiments, this mass shift will also change 
the number of formed J/ ip according to statistical hadroniza- 
tion. For T = 170 MeV, mass decrease of 100 MeV increases 
the yield by a factor of 2 since the factor is roughly charac- 
terized by e - Sm / T '. Furthermore, a hydrodynamic calculation 
shows the lifetime of QGP is ~ 4-5 fm/c at 200GeVM 
Au+Au collisions at RHIC 13211 . This will be much longer at 
LHC. If the width broaden to as large as such lifetime of the 
plasma, the mass shift might be detectable. 

In summary, we have given the first model-independent 
analysis of possible mass shift and broadening of width of 
J/tf) on the basis of lattice QCD inputs and QCD sum rule 
in the quenched approximation. Although the formalism, 
OPE up to dimension 4, is found to be applicable only to 
T ~ 1.05T C , we found that the change of gluon condensates 
in the deconfined phase causes notable reduction of the ratio 
of the moment r„|oPE, which results in mass shift and width 
broadening. The mass is found to be almost linear decrease 
with temperature if the width remains unchanged while the 
width linearly increases with temperature in the case of no 
mass shift. The maximum values are 5m ~ 200 MeV and 
r ~ 140 MeV at T = 1.05T C . Further details includingre- 
sults for different £ will be given in a future publication B23I1 . 
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